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1 INTRODUCTION 



In quantum cosmology instantons, solutions of the classical Einstein equations in Euclidean 
space, play an important role giving the main contributions to the path integral 0. Among 
them classical wormholes are of special interest, because they are connected with processes 
changing the topology of the models [0, |3|. We remind that classical wormholes are usually 
Riemannian metrics consisting of two large regions joined by a narrow throat (handle). 
They exist for special types of matter 0-0] and do not exist for pure gravity. In quantum 
cosmology it is generally assumed that on Planck scale processes with topology changes 
should take place. For this reason Hawking and Page H introduced the notion of quantum 
wormholes as a quantum extension of the classical wormhole paradigma. They proposed 
to regard quantum wormholes as solutions of the Wheeler-DeWitt (WDW) equation with 
the following boundary conditions: 

(i) the wave function is exponentially damped for large spatial geometry, 

(ii) the wave function is regular when the spatial geometry degenerates. 

The first condition expresses the fact that space-time should be Euclidean at spatial 
infinity. The second condition should reflect the fact that space-time is nonsingular when 
spatial geometry degenerates. For example, the wave function should not oscillate an 
infinite number of times. 

The given approach extends the number of objects which can be treated as wormholes 

II- 

We believe that for the description of quantum gravitational processes at high energies 
the multidimensional approach is more adequate. Modern theories of unified physical in- 
teractions use ideas of hidden (or extra) dimensions. In order to study different phenomena 
at early stages of the universe one should use these theories or at any rate models keeping 
their main characteristics. But more reliable conclusions may be done only on the basis of 
exact solutions which are usually obtained in rather simple cases. 

Therefore, at the beginning we consider a cosmological model with n (n > 1) Einstein 
spaces containing a massless minimally coupled scalar field and a cosmological constant 
A. The gauge covariant form of the WDW equation was proposed in [0. This model is 
integrable in the case with only one of the Einstein spaces being not Ricci-flat and vanishing 
cosmological constant. The general properties of this particular model were investigated in 
|| while classical as well as quantum wormhole solutions were found for different models 
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in || - fl2|| . (In [TJ| particular integrable cases with milticomponent perfect fluid were 
considered) . 

The present paper is devoted to the case of nonzero cosmo logical constant. For the 
model with one space of positive constant curvature in four dimensional space-time with 
cosmological constant and axionic matter (which is equivalent to a free minimally coupled 
scalar field, see, for example, |2|, [| and the paper by Brown et al. in [|]]) classical wormhole 
solutions were obtained in . In the case of four dimensional space-time with nontrivial 
topology R x S 1 x S 2 and non-zero cosmological constant this type of solution exists, too 

Here, we investigate the case with at least one of the spaces, say Mi being Ricci-flat. If 
all the other spaces M iy i = 2, . . . ,n are also Ricci-flat this model is fully integrable in the 
classical |12| as well as in the quantum cases @. For A < a family of quantum wormhole 
solutions with continuous and discrete spectra exist. Classical wormholes can be found in 
this case only in the presense of a scalar field. In the presence of a scalar field classical 
wormhole solutions exist also for another particular case with fine tuning of the parameters 
of the model, if A < and all Mj, i = 2, . . . , n are not Ricci-flat and have the same sign 
of the curvature. In this case only Mi has a dynamical behaviour and is considered as our 
external space. All the other internal spaces Mi,i = 2, . . . ,n are freezed with fixed scale 
factors a(o)j which are fine tuned to values determined by the cosmological constant. This 
type of solutions belongs to the class of models with spontaneous compactification. In the 
case of models without a cosmological constant and with only one non-Ricci-flat factor 



space solutions with spontaneous compactification were also found in fig] . 

We would like to note that solutions of the WDW equation in four dimensional models 
with A / and with a conformal scalar field were first obtained in [|T7| and |TB[ respec- 
tively (see also |0J). They include possibly the first quantum wormhole type solutions in 
four dimensions as well as DeWitt's solution for the Friedman universe with dust (1967). 
Vacuum quantum cosmological solutions in four dimensions may be found in f20| . The 
path integral approach to quantum cosmology [|IJ for models with cosmological constant 
in four and five dimensions with nontrivial topology was developed in [|2TJ . 

The paper is organized as follows. In section 2 the general description of the models 
considered is given. In section 3 classical and quantum wormholes are obtained for all 
spaces being Ricci-flat. In section 4 classical wormholes are considered in the model with 
spontaneous compactification of extra dimensions. Conclusions and an extensive list of 
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references complete the paper. 



2 GENERAL DESCRIPTION 
OF THE MODEL 

The metric of the model 

n 

g = ~exp[2-f(T)}dr ®dr + J2 exp[2P\r)]g {i) , (2.1) 

i=i 

is defined on the manifold 

M = 1R x Mi x . . . x M n , (2.2) 
where the manifold Mj with the metric is an Einstein space of dimension di, i.e. 

Rnnntb®] = (2-3) 

i — 1, . . . , n; n > 2. The total dimension of the space-time M is D = 1 + Yh=i c?j. 

Here we investigate the general model with cosmological constant A and a homogeneous 
minimally coupled field ip(t) with a potential U((p). 

The action of the model is adopted in the following form 

S = \ j d D xJ\g~\{R\g} - d M( pd N( pg MN - 2U{cp) - 2A} + S GH , (2.4) 

where R[g] is the scalar curvature of the metric g = gMN<lx M <8) dx N and S'gj/ is the 
standard Gibbons-Hawking boundary term |22j . The field equations, corresponding to the 
action fl2.4j) , for the cosmological metric (|2.1|) in the harmonic time gauge 7 = Y^l=\^S % 
are equivalent to the Lagrange equations, corresponding to the Lagrangian 

1 n 

L = g E (Gal?? + 2 ) - V, (2.5) 

with the energy constraint imposed 

1 n 

E=~I1 (GijFP 3 + v 2 ) + V = 0. (2.6) 

Here, the overdot denotes differentiation with respect to the harmonic time r. The com- 
ponents of the minisuperspace metric read 

Gij = dibij - didj (2.7) 
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and the potential is given by 

V=V(f3,v)=ex P (2jrd i p l 



i=i 



£>e- 2 ^ +U(<p)+A 



3=1 



(2.8) 



where 9i = \ l di. If the Mj are spaces of constant curvature, then may be normalized in 
such a way that 0j = kidi(di — 1), fcj = ±1, 0. We may also consider the generalization of 
the model with the potential fl2.8|) modified by the substitution 



(2.9) 



This gives us the possibility to investigate models with an arbitrary scalar field potential 
U(<p,P) = U(ip) as well as (for <p =const) models with an arbitrary potential U(ip,f3) = 
U(f3). Effective potentials of the form U(/3) may have their origin in an ideal fluid matter 
source. In special cases the general form U((f, (3) of the potential leads us to new integrable 
models. An example of this kind of potential will be presented. 
With the general potential U((p,f3) the equations of motion are 



k=l 

dU 



{^-1)9^ +d^e k e-^ k 



2d,U - 2d- A 







l,...,n, 



dU ( n N 
Cp +— exp (2^d# 



0. 



(2.10) 
(2.11) 



with the constraint 

At the quantum level the constraint (|2l|) is modified into the WDW equation (see fTj) 



2 I dp* dp d<J) 2 



(2.12) 



where \I/ = ^((3, f) is the wave function of the universe, V is the potential 

5 ij 1 

G lJ = — + 

<L 2 - D 



and 



(2.13) 



are the components of the matrix inverse to the matrix (Gy) ( |2.7| ). The minisuperspace 
metric G = G^dx % ® dx l (|2.7|) was diagonalized in []7], |23| 



G = -dz° ®dz + J2 dz i <g> dz { 



(2.14) 
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where 



j) -l 



n 

[di/H^] 1 ' 2 £ d0-(3% (2.15) 



i = 1, . . . , n — 1, where 



q=[(D-l)/{D-2)]W £ ? , = £^, (2-16) 

j=i 

The WDW equation (|2.12|) takes in variables ( |2.15|) the following form 

* = 0. (2.17) 



d d ^ d d d 2 n „, . 



dz° dz° ~[ dz % dz l dcp 2 



3 WORMHOLES FOR 
RICCI-FLAT SPACES 

In this chapter we consider the Ricci-flat case (9i = A^cf = 0, i = 1, . . . , n), with A / 
and U(ip) = 0. If the Mj are internal spaces they should be compact. This compactness 
is a necessary condition also for the Hartle-Hawking boundary condition (see below). The 
compactness of Ricci-flat spaces may be achieved by appropriate boundary conditions. The 
d-dimensional tore is the simplest example. 

3.1 Classical solutions 

In the considered case the Lagrangian ( |2.5|) may be written in the following form 

-i n+l n+1 

L = o E GuFF ~ Aexp(]T ujP 1 ), (3.1) 
z i,j=i i=i 

where [3 n+1 = ip, Ui = 2d{, i = 1, . . . , n, u n+ \ = and 

(Gu) = | ^ \ | (3.2) 
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(the matrix (Gij) is defined in eq. (|2.7|)). We consider the coordinates (z A ) = (z a ,z r ' 
/3 n+1 = <p) } where z a , a = 0, . . . , n — 1, are defined in ( |2.15| ). It is clear that 



A — 0, ... ,n, where 



n+l 

z A = Y, Vf?, (3.3) 

7=1 



1 

and the matrix (Vf) is defined in ( |2.15| ). This introduced matrix diagonalizes the minisu 



perspace metric 

n 

G u = £ VabVj A Vj B , (3.5) 

A,B=0 

I, J = 1, . . . , n + 1 and (r] AB ) = (r} AB ) = diag(-l, +1, . . . , +1). 

In the coordinates ( |3~3| ) the Lagrangian ( |3.1| ) reads (q is defined in ( |2.16|) ) 

1 n 

L =n Y, Vabz a z b -Aexp(2qz°). (3.6) 

Z A,B=0 

The Lagrange equations for the Lagrangian (|3.6| ) 

-z° + 2qA exp(2qz°) = 0, (3.7) 
z A = 0, A = l,...,n, (3.8) 



with the energy constraint 



1 n 

E =~Y1 Vabz A z b + Aexp(2qz°) =0 (3.9) 

1 A,B=0 



can be readily solved. First integrals of ( |3.8j ) are 

z A =p A , A=l,...,n, (3.10) 
where p A are arbitrary constants of integration. Then the constraint ( |3.9| ) may be rewritten 

- ^(z ) 2 + S + Ae 2qz ° = (3.11) 

with 

n 

2S=J2(P A ) 2 - (3-12) 

A=l 
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We obtain the following solution 

z A =p A r + q A , A = l,...,n, (3.13) 

where p A and q A are constants and 

2qz° = \n[£/{Asmh 2 (qV2£(r - r ))}], £ ^ 0, A > 0, (3.14) 
= ln[l/{2g 2 A(r-r ) 2 }], £ = 0, A > 0, (3.15) 
= \a[-£/{Acosh 2 {qV2£{r -T ))}\, £>0,A<0, (3.16) 

Here t$ is an arbitrary constant. 

3.2 Kasner-like parametrization 



First we consider the case £ > 0. In this case the relations (3.14) and ( |3.16| ) may be written 
in the following form 



2qz° = \n[£/{\A\f 2 (qV2£(r - r ))}], (3.17) 
where 5 = A/|A| = ±1 and 

fs(x) = i(e x — 5e~ x ) = sinhx, 5 = +1, 

= coshx, 5 = -I. (3.18) 

We introduce a new time variable by the relation 

t= r iD expfaV2£(r-r )) + ^ p lg) 
V<5 exp(gv2£(r — r )) — v 5 

where 

T= [( J D-2)/2|A|(D-l)] 1 / 2 . (3.20) 
It is not difficult to verify that the following relations take place 

sinh(tVd/T)/Vd = l/f s {q^2£{r - r )), (3.21) 

tanh(tV5/2T)/V5 = exp{-qV2£{r - r )), (3.22) 

dt = -TqV2£dr/f s (qy/2£(T-T )). (3.23) 

Now, we introduce the following dimensionless parameters 

n 

oi 1 = ~ E Vlp A /qV2£, (3.24) 

A=l 
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where (Vj) = (V/ 4 )" 1 (see Q). It is clear that 

(n) - 



v: o 



1 



(3.25) 



where (V a l ) = (V?) 1 . The relation (|3.24j) is equivalent to the following relations 

n-l _ 

-EKVV^, i = l,...,n, (3.26) 



a 1 



a 



a=l 
n+l 



-p n /qV2£ 



(3.27) 



It follows from eq. (|3.5|) , that 



n+l n 

V{ = E E G IJ Vfr]BA 

,J=1 B=0 



where 



= {Guy 1 = 

From ( |2.13D and ( |3.28| ) we have 

V* = -G^ Uj /2q = (q(D - 2))" 1 



G ij 
1 



V " +1 = 0. 



Using the relations ( pT3|) , (|3~T7D , O)-(^), ( gjg) and (|3T30|) we get the 
expressions for the solution of field equations 



g = -dt (g)dt + Y^ a 2 i(i)9{ 



i=l 



Oi(t) = exp(/?(t)) = A i [smh(rt/T)/r] a [tanh(rt/2T)/rf 1 
exp(<^)) = exp(^ +1 (t)) = A n+1 [tanh(rt/2T)/r]' 5 " +1 , 



where t > 0, r = y/6 = y A/|A| = V±T, a — (D — 1) 1 , 7^ are constants, i 
and the parameters a 1 satisfy the relations 

1 n+l n 

- e «/« / = E^« l = °> 
1 1=1 i=i 

n+ 1 n 

E G 7J d 7 d J = E d ^f + («" +1 ) 2 = ( D - 2 )/( D - !)• 
i,j=i i=i 



(3.28) 
(3.29) 

(3.30) 
folowing 

(3.31) 

(3.32) 
(3.33) 



l,...,n, 



(3.34) 
(3.35) 
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The first relation ( |3.34| ) can be easily proved, using the definition ( |3.24| ) and the following 
identity (we remind that due to fl2.15| ) and ( ^.3| ) uj = 2qVj) 

n+l n+l 

J2 UiV{ = 2g ]T VfV{ = 2q5° A = (3.36) 
1=1 i=i 

for A > 0. The relation ( ggg) follows immediately from ( p32p , ( ggg) and 

(PI). 

Now we consider the case £ = 0. In this case for A > there exist also an exceptional 



solution with the following scale factors in (p. 31 



ai{t) = Aiexp(±at/T). (3.37) 

i = 1, . . . ,n, and ip(t) =const. (This solution can be readily obtained using the formulas 
( gl3| ) and (gl5|) .) 

It is interesting to note that for A > the solution ( [3.371 ) with the sign "+" is an 
attractor for the solutions (13.321), i.e. 



di(t) ~ Ai exp(crt/T), i = l,...,n, (3.38) 

and ip(t) ~const for t — > +oo. The relation ( |3.38|) is the isotropization condition. We 
note that the solution (|3T3TD - (|3T35|) with a n+1 = was considered previously in |p6| . The 
special case of this solution with n = 2 was considered earlier in J?7] . 
The volume scale factor corresponding to (|3.32j) has the form 



n n 

v = H at = ([] Af) sinh(rt/T)/r (3.39) 
i=i i=i 

It oscillates for negative value of cosmological constant and exponentially increases as 
t — > +oo for positive value. For positive Ai and E the following identity takes place 

n 

i=l 



l[At = ^£/\A\. (3.40) 
i=i 

For small time values we have the following asymptotical relations 

Oi(t) ~ Ci t ai , exp^(t) ~ c„ +1 t Q " +1 (3.41) 
as t — > 0, i = 1, . . . , n, where 

on = ct + a, a n+1 = a n+1 , (3.42) 



10 



are Kasner-like parameters, satisfying (see ( |3.34| ), ( |3.35| )) the relations 



ai) 2 + a 2 n+1 = 1. 



(3.43) 



i=l 



i=l 



The behaviour of the scale factors near the singularity coincides with that for the case 



A = |§, gg]. For a n+1 = see also p 



We note also that in terms of arparameters the solution ( |3.31| ) - (|3.35|) reads 

ai(t) = A t [smh(rt/2T)/r}^[cosh(rt/2T)} 2a ~ a \ (3.44) 
exp(<p(t)) = A n+ i[tanh(rt/2T)/r] a " +1 . (3.45) 



Let us apply these solutions to the case of two spaces (n = 2). From ( |3.43| ) and ( |3.44| ) we 
find for the non-exceptional solutions 



a 2 



Ai[sinh(rt/2T)/r] a i [cosh(rt/2T)] c 
A 2 [sinh (rt / 2T) /r) a * [cosh (rt /2T) ] c 



where 



and 



(X i 



dt ± VR 
d 1 (d 1 + d 2 ) 



a .-> 



d 2 ±^R 
d 2 (di + d 2 ) 



(3.46) 
(3.47) 

(3.48) 
(3.49) 



R = cM 2 [K + d 2 )(l - aj) - 1]. 

Graphically these solutions are presented in Figs. 1-2. 

In the Euclidean case after the Wick rotation (t — > —it) we get the following instanton 
solutions 



g = dt <g> dt + XX ? (*)0(O> 
i=i 

Oi(t) = i i [sinh(ts/2T)/s] Q! [cosh(ts/2T)] 2<T - Q % 
exp((p(t)) = i„+i[tanh(ts/2T)/s] a " +1 . 



(3.50) 

(3.51) 
(3.52) 



where T is defined by ( |3.2(J| ), s = J— A/|A| and the parameters ai satisfy the relations 
( p.43| ). For A < we have the special solution (£ = 0) 



di(t) = Aiexp(±at/T). 



(3.53) 



We note that for the Euclidean case the scale factors may be obtained from the correspond- 
ing Lorentzian ones by the substitution A i— > —A. For n = 3, d\ = d 2 = d% = I, A > the 
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special case of the solution ( |3.50| ) - ( |3.52| ) was considered in [p0|] . We note that for A < 
there are wormhole-like sections of the total metric ( |3.50| ). This takes place, for example, 
if n — 2, a| < 1 — d2 , 1 < d\ < d 2 , (see Fig. 2). In this case the scalar field is real in 
Euclidean region. 

Now, we consider the solutions of the field equations with complex scalar field and real 
metric. In this case £,p 1 , . . . ,p n_1 are real and hence (see ( |3.12| )) p n is either real or pure 
imaginary. The case of real p n was considered above. 

For pure imaginary p n we have three subcases: a) £ > 0, b) £ = 0, c) £ < 0. In the first 
case a) £ > after the reparametrization ( p,19|) , (|3.20|) we get the solutions ( |3.32|) -( |335"D 
with an imaginary value of a n+1 . The cases b) and c): £ < take place only for A > 0. 



Let us consider the case c) £ < 0. Here, we have (see ( |3.26| ), ( |3.27| )) imaginary 



a 



i a k , k = 1, . . .,n, at'" - = cr n+1 . 
The solution may be obtained from ( ^.32[ )-( p75o] ) substituting ( |3.54| ) and t/T 

n 

g = -dt(g)dt + J2 a2 i > 

a,i(t) 

<p(t) 



(3.54) 



t/T 



2 ' 



i=l 



i 4 [cosh(t/T)] CT [/(t/2T)]^, 
c + 2ia n+ i arctan e~ t//T , 



(3.55) 

(3.56) 
(3.57) 



where c, A, ^ are constants, i — 1, . . . , n, a — (D — 1) 1 , T is defined in ( |3.20|) , A > 
and the real parameters o - / satisfy the relations 



n 

- d i°1 + °"n+l 



0. 



(D-2)/(D 



i=l 



Here 



,7T, 



/(x) = [tanh(a; + z^)] 4 = exp(-2 arctan e 2x ) 



(3.58) 
(3.59) 

(3.60) 



is smooth monotonically increasing function bounded by its asymptotics: 
e _7r < f(x) < 1; f(x) — > 1 as x — > +oo and /(a;) — > e _7r as x — >• — oo (see Fig. 3). 
The solution ( 3.55|) - (|3.59|) may be also obtained from formulas ( |3.13|) , ( p.!4j ). The relation 
between the harmonic and the proper times ( |3.21| ) is modified for our case £ < 



cosh(t/T) = l/sm(qJ2\£\(r-r )). 



(3.6i; 
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For the volume scale factor we have 



n 

v = IT " 

i=l 



di 



l[A?)cosh(t/T). 



(3.62) 



i=i 



The scalar field varies ip(t) varies from c + ma n+ i to c as t varies from — oo to +00. The 
solution ( |3.55| )- (|3.59|) is non-singular for t 6 (—00, +00). Any scale factor a^t) has a 
minimum for some to« an d 

Oi(t) ~ Afexp(a\t\/T), 



(3.63) 



for t — > ±cx). 

The Lorentzian solutions considered above have also Euclidean analogues for 
A < 



g = dt £g> dt + crf(t)g( 



0' 



a,(t) = i l [cosh(t/T)] CT [/(t/2T)]^ 
y?(t) = c + 2ia n+ i arctan e _< ^ T , 



(3.64) 

(3.65) 
(3.66) 



with the parameters 07 satisfying the relations ( |3.58|) - (|3.59|) . This solution may be inter- 
preted as classical Euclidean wormhole solution. An interesting special case of solution 
( |3.64| )-( pi3(f ) occurs for 0{ = 0, i = 1, . . . , n, (this corresponds to p l = 0) 



Oi(t) = AilcosHt/T)]*, 

(p(t) = c ± 2iq~ l arctan e _ ^ T . 



(3.67) 
(3.68) 



All scale factors ( |3.67| ) have a minimum at t — and are symmetric with respect to time 
inversion: t 1— > — t. We want to stress here that wormhole soluions take place only in the 
presence of an imaginary scalar field in the Euclidean region. Analytic continuation of the 
solutions (|3.67|) , (|3.68|) into the Lorentzian region leads to real geometry and real scalar 
field there. 



3.3 Quantum wormholes 

The model introduced above leads to the WDW equation (|2.17|) 



2H^ 



d d n d d 

1=1 



^ = 0. 



(3.69) 



13 



We are seeking the solution of ( |3.69| ) in the form 

= exp{ipz)$(z°), (3.70) 

where p = (p 1 , . . . ,p n ) is a constant vector (generally from C n ), z = (z 1 , . . . , z n ~ x , z n = if), 
pz = YHl=\V% z% an d Pi = Yl 1 j=iVijP' = P The substitution of (|3.70| ) into (|3.69|) gives 

, 1, d 



l o , ()30 ) 2 + V (/)]$ = ^, (3.71) 

where £ = \pp and Vo(z°) = — Ae 2qz ° . The potential Vo(-2°) is plotted on fig. 4 and fig. 
5 for A > and A < respectively. The clasiically allowed (Lorentzian) and forbidden 
(Euclidean) regions are shown there with respect to the energy levels £. Solving ( |3.71| ), we 

get 

= 5v2£/ g (V=2Ag-V* ), (3.72) 
where i\pl£ jq = i\p\/q, and B = J, K are modified Bessel functions. We note, that 



n 



v = expqz = Y[a? (3.73) 

i=i 

is proportional to the spatial volume of the universe. 

The general solution of Eq. Q3.69D has the following form 



¥(*)= £ J d n p C B (p)e i ^B im/q (\f^2Aq- 1 e qz °), (3.74) 



B=I,K 



where functions Cb (B = I,K) belong to an appropriate class. Similar solutions were 
found for the two-component model (n = 2) and A > in [ fn|j . 

The solutions (|3.70|) are the eigenstates of the quantum-mechanical operators fl z i = 



— (i/N)d/dz l ,i = 1, . . . ,n with the eigenvalues (1/N)pi where N = 1 for the Lorentzian 
space-time region and N = i for the Euclidean one. 

Due to the well known time problem in quantum cosmology the WDW equation is not 
really the Schrodinger equation. There is no generally accepted procedure to overcome 
this problem but for our particular model we can introduce some time coordinate into the 
quantum equations in analogy to ||. 

We split the WDW operator H ( p. 69 ) into two parts 



H = -H + H u (3.75) 
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where 



and 



1 d 2 



Ae 2 < 



1 " d 2 



(3.76) 
(3.77) 



Then the WDW equation ( |3.69| ) becomes 

H ^ = Hi®. (3.78) 

Applying Hi to the wave function ( p.70| ) one gets 

= EV. (3.79) 

Now, we take £ to be real. Then, equation ( |3.79| ) shows that £ can be treated as the 
energy of the subsystem Hi and equation Q3.79D becomes the Schrodinger equation. From 
this point of view ^ gives the stationary states of the subsystem described by the wave 
equation (in Lorentzian region) 

i— = Hi^>, (3.80) 

where 

ij> = e -^ r ^. (3.81) 
It can be easily seen that the wave equation 



i— = H Q m 

OT 



(3.82) 



is reduced to equation (|3.71| ). 

In the semiclassical limit for the wave function ( |3.81 ) equations ( |3.80| ) and ( 3.82| ) are 
reduced to the classical equations ( |3.1U| ), ( P-llp . Indeed, the wave function ( |3.81| ) can be 
rewritten in the form 

§ = e~ i£r e iSl <5>{z°) (3.83) 

with Si = J27=iPi z ' ''■ m the semiclassical limit the wave function takes the form 

<$>(z°) = C{z°)e iSo (3.84) 

with C(z°) being a slowly varying function and So(z°) being a rapidly varying phase. Time 
is defined in the semiclassical limit as an affine parameter along integral curves 



d _™ d(s + Si) d 

Or dzi dz l 



(3.85) 
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where Zq 



-Z , Z{ 



1, ... ,n, As result we find the equations z % = p % 



1 n, 



and these coincide with the classical equations (|3.10|). For this reason we used the same 



notation for the constants of integration in (|3.10|) and for the momenta in the wave function 
( p.70| ). The velocity along z° is found to be i = — i° = dS /dz°. Using this relation and 



putting the wave function (|3.83|) , ( |3.84|) into equation ( |3.82| ) we reproduce the classical 



equation (3.11 



As shown above the parameter S can be interpreted as energy. So we may treat the 
state £ = as the ground state of the system. The demand of reality of the geometry 
leads to real momenta p* (i = 1, . . . , n — 1) in the Lorentzian region. The scalar field can 
be real or imaginary there. In the ground state we put all momenta p l (i = 1, . . . ,n) equal 
to zero and the ground state wave function reads 



*o = B ( g - 1 V=2Ae^°) . 



(3.86) 



It is interesting to note that is invariant with respect to the rotation group O(n) in the 
space of vectors z = (z 1 , . . . , z n ). 

In eq. ( 3.86 ) Bq denotes the Bessel functions of order zero. A particular solution may be 
specified by boundary conditions. For example, quantum wormhole boundary conditions 
were presented in the Introduction. Among the different types of boundary conditions 
for wave functions describing the universe the most popular is the Hartle-Hawking (HH) 
boundary condition M. According to the HH proposal the ground-state wave function of 
the universe is given by a path integral over all compact Euclidean geometries and 
the regular matter fields: 



{HH) 



d[g]d[ip]e 



-Is 



(3.87) 



where Ie is the Euclidean action. For our model the Euclidean action in harmonic time 
gauge and in z-coordinates reads 



I E = - [ T dr \-(z ) 2 + + 2Ae 29Z ° 

2 Jr * I i=i 



It) 
2v 



(3.1 



where v denotes the spacial volume ( |3.73| ) up to a numerical factor. The upper limit r 
corresponds to the boundary of the D-dimensional manifold, where the z % {% = 0, . . . , n) 
have values indicated by the arguments of . The lower limit r* corresponds to the 

point where the D-dimensional manifold closes in a smooth way. The origin of the second 
term in ( |3.88j ) was explained in detail, e.g. by Louko In the semiclassical limit the 
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wave function is given by 

tt^we-'*, (3.89) 

where Jg should be calculated on the classical Euclidean solutions with boundary condi- 
tions defined by the concrete scheme of geometry closing at r = r*. 

Now, we find the relationship between the HH wave function (|3.87|) , ( |3.89| ) and our 
ground-state wave functions (|3.86j ). Let us first consider the case of a negative cosmological 
constant A < 0. Then, we have the classical Euclidean equations 



i u = ±J2\A\e qz . (3.90) 



The spacial volume v may be presented in the form 

-l 



v = e qz =- [q^2\A\ t) , -oo < r < 0. (3.91) 

Formula ( |3.91|) shows that the geometry closes at the harmonic time r — > — oo. It is easy 
to see from ( |3.91| ) that the second term in ( |3.88| ) contributes nothing to Ie- So, on these 
classical solutions the Euclidean action reads 



q* t q 

and we get the semiclassical HH wave function ( |3.89| ) 



1 1 -a/2 A o 
J e = Z2Z = -^—-e** (3.92) 



- exp \l^ e ^ ' j = exp ( ^ g <#) ■ ( 3 - 93 ) 

Eqn. Q3.93|) shows that (for the class of real Euclidean geometries) 

® { o HH) ^ + oo, z°^+oo. (3.94) 



This conditions provides us with the possibility to chose a solution of equation ( 3.69 ) 
corresponding to the HH ground state: 



Q 

The vacuum solution ( |3.95| ) has the asymptotic form 



(3.95) 



*i HH U exp [ V^e^° 1 , z ^ +oo, (3.96) 
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which coincides with ( |3.93| ). 

A similar procedure can be performed for a positive cosmological constant A > (see, 
e.g. f|). In this case the classical Euclidean equation is 



+ 2q 2 Kv 2 = 



(3.97) 



and gives an imaginary geometry. This reflects the fact that the geometry should be 
purely Lorentzian in the case A > for £ > 0. The action I|? ( |3.88| ) is indefinite in this 
case. We can avoid this problem if we perform the analytic continuation v — > iv. After this 
continuation the action ( 3.88Q is formally reduced to the action in the case A < 0. Again, it 

( HH) 

leads to the following asymptotic for the HH wave function: \l/g ^ + oo, v — > +oo, which 
shows that the wave function I [(y/2A/q)v] is a solution of eq. (|3.69| ) corresponding to the 
HH ground state for the class of Euclidean solutions considered here. Thus, after analytic 
continuation to real values of v the vacuum state corresponding to the HH boundary 
condition is 



(HH) 



This solution has the asymptotic form 



Jn 



(HH) 



2A 



cos 



COS 



d, 



+ 00. 



(3.98) 



(3.99) 



i=i 



For the Bianchi I universe (n = 3, d\ = d 2 = d 3 = 1) eq. (|3.99|) is reduced to 



^ { HH) ~cos (2JA/3a 1 a 2 a 3 



+oo. 



(3.100) 



Similar results for the Bianchi I universe were obtained earlier in papers by Laflamme and 
Louko 0]. 

Now, let us turn to quantum wormholes. We restrict our consideration to real values 
of Pi. This corresponds to real geometries in the Lorentzian region. In this case we have 
£ > 0. 

If A > the wave function \& ( 3.70| ) is not exponentially damped when v — > oo, i.e. the 
condition (i) for quantum wormholes (see the Introduction) is not satisfied. It oscillates 
and may be interpreted as corresponding to the classical Lorentzian solution. 

For A < 0, the wave function ( |3.70| ) is exponentially damped for large v only, when 
B = K in ( |3.72| ). But in this case the function $ oscillates an infinite number of times, 
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when v — > 0. So, the condition (ii) is not satisfied. The wave function describes the 
transition between Lorentzian and Euclidean regions. 
The functions 

tfiKz) = e^K^iV^Xq- 1 ^* ), (3.101) 

may be used for constructing quantum wormhole solutions. Like in []S], |TU|, we consider 
the superpositions of singular solutions 

$>x,n( z ) = ~ / dkS> qkR (z)e- ik \ (3.102) 

71 J-oo 

where A G 1R, n is a unit vector (n 2 = 1) and the quantum number k is connected with 
the quantum number £ = | | p \ 2 by the formula IE = q 2 k 2 . The calculation gives 

#A,tfO) = exp[-^-- — e 9 *° cosh(A - gzn)]. (3.103) 

It is not difficult to verify that the formula ( |3.103| ) leads to solutions of the WDW equation 
( 3-69| ), satisfying the quantum wormholes boundary conditions. 



We also note that the functions 

( r 0\2 , / 1\2 

* m>H = H m (x°)H m (x l ) eM~ } a ^ ], (3-104) 

where 

x° = (2/g) 1/2 (-2A) 1/4 exp(g//2) cosh(^5n), (3.105) 



x 1 = (2/g) 1/2 (-2A) 1/4 exp(g//2) sinh(^£n), (3.106) 

m — 0,1, ... , are also solutions of the WDW equation with the quantum wormhole bound- 
ary conditions. Solutions of such type were previously considered in 0, |TU). They are 
called discrete spectrum quantum wormholes. 

It is clear from the equation (|3.71|) and fig. 4 that in the case A < a Lorentzian region 
exists as well as an Euclidean one for £ > 0. In the case A > only the Lorentzian region 
occurs for £ > and for £ < both of these regions exist (see fig. 5). The condition £ < 
leads for pure gravity to a complex geometry in the Lorentzian region. We can avoid this 
problem by the help of a free scalar field, because in this case 2£ = Y^Zi pf + Pn and we 
can achieve £ < for real pi(i = 1, ... ,n — 1) and imaginary p n , i.e. for an imaginary 
scalar field in the Lorentzian region. The wave functions ( |3.70|) , ( |3.72| ) with A > and 
£ < describe the transitions between Euclidean and Lorentzian regions, i.e. tunneling 
universes. 
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4 CLASSICAL WORMHOLES, 

FINE TUNING OF PARAMETERS 



Classical wormhole solutions exist in our model also for another interesting case. This 
is the case with spontaneous compactification of the internal dimensions. Let the factor 
space Mi be our dynamical external space. All the other factor spaces Mj(z = 2, . . . , n) 
are considered as internal and static. They should be compact and the internal dimensions 
have the size of order of Planck's length Lpi ~ 10~ 33 cm. The scale factors of the internal 
factor spaces should be constant: Oj = = O(o)i(i — 2, . . . ,n). It is not difficult to show 
that in the case of fine tuning of the parameters due to 

ft 2A 



d i a %)i D ~ 2 



Co, i = 2,...,n, (4.1) 



all dynamical equations (|2.10 ) are reduced to one for the scale factor ai = and this 
equation reads 



The constraint (??) has form 



(4.2) 



di (di-l)$ 12 =^-eie 2 ^-^\ 2 ^ k +e 2d ^e^ dkPk [2U+ 

+C (di - 1)] . (4.3) 

From the equations (|4. 1|) it follows that all internal spaces should be non-Ricci-flat and 
signft = signA, (i = 2, . . . , n). We remind here that overdot denotes differentiation with 
respect to the harmonic time r 0. The minimally coupled scalar field has the specific 
potential U(/3,ip) = U(<p)exp [— 2 Yh=i di{3 1 } ■ For this potential it is easy to get the first 
integral of equation ( |2.11|) 

if 2 + 2U(ip) = v 2 = const. (4.4) 



This gives 

for U (<£>) = and 



<f — ut + const (4.5) 
if = (fo cos m(r — r ) (4.6) 
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for Ufa) 



where v = nupo here. 



Let us investigate the model where our external space Mi is Ricci-flat, i.e. 9i = 0. 



Then we can rewrite the equation (|4.3| ) as follows 

1 ) 2 = 9 2 + ~Ae 2d ^ , 



where the constants are 



and 



v 2 



A 



&i{di 
2A 



tt5 



l (0)fc- 



(4.7) 



(4.9) 



^1 (Efc=l dk fc=2 

It is clear from equation (^4.7|) that the dynamical behavior of the scale factor ai depends 
on the signs of v 2 and A . If A > and v 2 > then ai expands from zero to infinity. 
For A > and v 2 < 0, ai has the turning point at some minimum and this case may be 
realized for an imaginary scalar field in the Lorentzian region. For a real scalar field in the 
Lorentzian region ( i.e. v 2 > ) and A < the scale factor ai expands from zero to its 
maximum and after the turning point shrinks again to zero. For the latter case the solution 
has a continuation into the Euclidean region with the topology of a wormhole, that means, 
two asymptotic regions which are connected with each other through a throat. 

Let us investigate the case with A < in more details. As signA = sign#j (i — 2, . . . , n) 
then for A < the curvatures 0, < also. (As a special case the internal spaces Mj 
(i — 2, . . . ,n) may be compact spaces of constant negative curvature [[25j.) The solution 
of equation ( [4.7| ) (the Lorentzian region) has the form 



aAr) 



n ~ 1 l/2di 

v /|A| [cosh di vt 



-l/dl 



-oo < r < +oo. 



(4.10) 



The synchronous time t and the harmonic time r are connected by the differential 
equation 

e^dr = dt, (4.11) 

where 

n 

7(r)=5>0*- (4.12) 



i=l 



It is not difficult to get the connection 

cosh(<iiZ>T) = 



cos ( \l Codi t + const J 



i -i 



(4.13) 
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With the help of this connection we obtain the expression for the scale factor eti with 
respect to the synchronous time 

I l/di 



ai(t)= ^/|A| 



l/2di 



sin \/C di t 



0<t< 



71 



VCodi 
0. For i 



(4.14) 
we have 



where the constant in ( |4.13|) was fixed by condition a(t = 0) 
<2i ~ t 1 ' 1 * 1 . Thus the external space Mi has the behavior of a FRW-universe filled with 
radiation for d\ — 2 and with ultrastiff matter for d\ — 3. The Lorentzian metric ( |2.1|) in 
synchronous time gauge reads as 



g = -dt®dt + a\(t)g(X) + a to)i9(i) 



(4.15) 



8=2 



with ai given by ([4.14|) . 

Our next step is to get the wormhole-type solution for Mi in the Euclidean region. The 
transition into the Euclidean space is performed by the Wick rotation t — > —it. The exact 
form of the transformation from the Lorentzian time to the Euclidean "time" can 
be obtained demanding the existence of wormholes being symmetric with respect to the 
throat (see Zhuk in [[|). In what follows, for the expression ( |4.14j ) we should perform the 



analytic continuation 



ai(f)= *7|A| 



2 V / CWI 

l/2di 



— He- This gives us 

l/di 



cosh w Cod\t 



-oo < t < 



-oo. 



(4.16) 



Of course, this formula can be obtained also as a solution to the Euclidean analog of the 
equations (|4.3|) , (|4.7| ) for an imaginary scalar field in the Euclidean region. The metric of 
the Euclidean region is given by 



g = dt®dt + a\{t)g(X) +J2 a lo)i9(i), 



(4.17) 



i=2 



where a x (t) is described by fl4.16|) . Thus, in Euclidean space we have two asymptotic 

r ~ -I l/2di 

regions t — > ±oo connected through a throat of the size v /|A| and this object is a 
wormhole by definition. 

It is clear that the Lorentzian solution ( [4. 14j ) and its Euclidean analog ( |4.16| ) take place 
only in the presence of a real scalar field in the Lorentzian region (i.e v 2 > 0) or equivalently 
an imaginary scalar field in the Euclidean region. 
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5 CONCLUSIONS 

In this paper we investigated multidimensional cosmological models with n(n > 1) Einstein 
spaces Mi in the presence of the cosmological constant A and a homogeneous minimally 
coupled scalar field ip(t) as a matter source. The problem was to find classical and quantum 
wormhole solutions. Classical wormholes are solutions of the classical Einstein equations 
describing Riemannian metrics with two large regions joined by a throat. Quantum worm- 
holes are solutions of the Wheeler-DeWitt (WDW) equation with the proper boundary 
conditions proposed by Hawking and Page 0. 

The model was investigated where one of the factor spaces, say Mi, is Ricci-flat. In the 
case when all other factor spaces Mi, % = 2, . . . , n, are Ricci-flat too, the classical Einstein 
as well as the WDW equations are integrable. For a negative cosmological constant A < 
quantum wormhole solutions were constructed. These solutions exist for pure gravity 
as well as for the model with a free minimally coupled scalar field. Classical wormhole 
solutions exist in the Euclidean region for A < in the presence of an imaginary as well 
as real scalar field. 

Classical wormhole solutions were also obtained in models with spontaneous compacti- 
fication. In this case the Ricci-flat factor space Mi was considered as our external dynam- 
ical space. All other factor spaces Mj, i = 2, . . . , n are static with constant scale factors 
a {o)i — const and all of them are fine tuned to each other and to the cosmological constant: 

Oj — §k — _2A_ „• u _ o „ 

7777T 2 r>—i > > > • • • > 

a ^ a (0)i aka (o)k u z 

As in the previous model, wormhole solutions exist for a negative cosmological constant 
A < 0. But there are important differences. Firstly, all inner spaces Mj,i = 2, . . . ,n, are 
non-Ricci-flat and have negative curvature. Secondly, the wormhole solution for the later 
case exists only in the presence of an imaginary sclar field in the Euclidean region. Thirdly, 
it seems hardly to be possible in the case 82, ■ ■ ■ , 9 n 7^ to integrate the Einstein equations 
as well as the WDW equation without the demand of spontaneous compactification with 
fine tuning. 

In models with one scale factor having a turning point (at the minimum) the production 
of the Lorentzian space-time is treated as a quantum tunneling process |32| ("birth from 
nothing"). The universe appears spontaneously going through the potential barrier with 
size equal to the size of the Lorentzian universe at the turning point. In our case of 
multidimensional models this kind of interpretation becomes more complicated. It follows 



23 



from ( |3.56| ) that the factor spaces Mi in general reach their minimum expansion positions 
at different times. The "birth from nothing" for each factor space takes place at a different 
value of time. If the difference between these events goes to infinity the extra dimensions 
are in the classically forbidden region forever. This interpretation is in the spirit of the 
Rubakov-Shaposhnikov idea |33| stating that extra dimensions are unobservable because 
they are hidden from us by a potential barier. 
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Figures: 



Fig. 1 The behaviour of the scale factors af(t),af(t) (see ( |3.46|) , (|3.47|) ) for A < 0, n = 
2, 1 < d\ < d,2 and different values of the integration constant a^. The left figures corre- 
spond to the sign "+" in ( ggp , ( ^47|) . 

Fig. 1.1 < a\ < 1 - i 



Fig. 1.2 a§ = 1 - £ 



Fi §- 1.3 1 - i < <4 < 1 - i 
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Fig. 1.4 «§ = 1 - i 



Fig- 1-5 1 - i < o§ < 1 - di+d2 



Fig. 2 The behaviour of the scale factors af(t),af(t) (see ( ET35D , (|3T47| )) for A > 0, n = 
2, 1 < <ii < d 2 and different values of the integration constant a 3 . The left figures corre- 
spond to the sign "+" in ( [Up , ( g747|) . 



Fig. 2.1 < a§ < 1 - £ 



Fig. 2.2 a\ = 1 - £ 
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Fig. 2.4 



a: 



Fig. 3. Representation of the function f(x) = exp(— 2 arctane 2x ). 



Fig. 4. The potential V (z°) for A > (solid line) and the energy levels E (dashed lines). 
For S > the Lorentzian region exists, only. For S < both regions, the Lorentzian as 
well as the Euclidean one exist. In this case quantum transitions with topology changes 
take place (tunneling universe). 
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Fig. 5. The potential Vq(z°) for A < (solid line) and the energy levels S (dashed lines). 
For S < the Euclidean region exists, only. For S > both regions, the Lorentzian as well 
as the Euclidean one exist. In this case quantum transitions with topology changes take 
place (quantum wormholes). 
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